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BACKGROUND OF THE INVENTION 

The present invention relates to floating-point processors, and more 
particularly to floating-point processors with an operating mode that has improved 
accuracy and high performance. 
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In digital processing systems, numerical data is typically expressed using 
integer or floating-point representation. Floating-point representation is preferred in 
many applications because of its ability to express a wide range of values and its ease of 
manipulation for some specified operations. A floating-point representation typically 
5 includes three components: a sign bit (sign), a mantissa (mant) that is sometimes referred 
to as a significand, and an exponent (exp). The represented floating-point number can be 
expressed as (-l) slgn »mant*2 exp . Floating-point representations are also defined by "EEEE 
Standard for Binary Floating-Point Arithmetic," which is referred to herein as the EEEE- 
754 standard (or simply, the IEEE standard) and incorporated herein by reference in its 
1 0 entirety for all purposes. 

The IEEE standard defines representations for "normalized" and 
"denormalized" float-point numbers. A normalized floating-point number is 
characterized by a mantissa having a one to the left of the binary point and a Lxxx— xx 
format, where each "x" represents one bit that is either a one or a zero. A denormalized 
Jj 1 5 floating-point number is characterized by a mantissa having a zero to the left of the binary 

I V point and a format of O.xxx—xx. Floating-point numbers greater than or equal to a 

i 1 

rh positive minimum representable normalized number (i.e., y > +a m j n ) and less than or equal 

JU = to a negative minimum representable normalized number (i.e., y < -a m i n ) are represented 

:i\ using normalized numbers. Floating-point numbers in the range between the negative 

fij 20 and positive minimum normalized numbers other than zero (i.e., -amin < y < +&mm, and y * 
3? ±0) may be represented using denormalized numbers. Zero is represented by a mantissa 

having a value of zero and an exponent also having a value of zero. 

Many operations can be performed on floating-point numbers, including 
arithmetic operations such as addition, subtraction, multiplication, division, and 
25 reciprocation. For arithmetic operations, the IEEE standard provides guidelines to be 

followed to generate a unique answer for each floating-point operation. In particular, the 
IEEE standard describes the processing to be performed on the result from a particular 
operation (e.g., add, multiply), the precision of the resultant output, and the data format to 
be used. For example, the IEEE standard defines several rounding modes available for 
30 the results from add and multiply operations, and the bit position at which the rounding is 
to be performed. The requirements ensure identical results from different 
implementations of IEEE-compliant floating-point processors. 
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The rounding modes defined by the IEEE standard provide improved 
accuracy for some operations, but are complicated to implement and also increase the 
processing time for an arithmetic operation. To obtain an output that fulfills IEEE 
rounding requirements, post-processing of a preliminary result from an arithmetic 
5 operation is typically performed. The post-processing includes possible denormalization 
and rounding of the preliminary result in accordance with one of the rounding modes 
defined by the IEEE standard. Denormalization is performed on a number having an 
absolute value less than +a min (i.e., -a min < y < +a min ) to place it in a proper format such 
that rounding can be performed at the bit location specified by the IEEE standard. The 
10 post-processing (or more specifically, the denormalization and rounding) typically leads 
to increased circuit complexity and increases processing time. The IEEE rounding modes 
are thus generally implemented in applications requiring high accuracy. 
]] To reduce hardware complexity and improve processing time, many 

i floating-point processors implement an operating mode in which numbers (e.g., 

?l 15 preliminary results) within the range of negative and positive minimum normalized 
j numbers (i.e., -a m j n < y < +a m i n ) are set or flushed to zero. The "flush-to-zero" mode is 

simple to implement and only marginally increases the processing time. However, the 
flush-to-zero mode suffers a loss in accuracy since the mantissa is flushed to zero. 

For many applications such as embedded processors, reduced cost and 
20 improved processing time are desirable. For these applications, an operating mode that is 
simple to implement, marginally increases the processing time (if at all), and has 
improved accuracy over the flush-to-zero mode is highly desirable. 



SUMMARY OF THE INVENTION 

25 The invention provides an operating mode (referred to herein as a "flush- 

to-nearest" mode) that has improved accuracy over a conventional "flush-to-zero" mode. 
In one embodiment, the flush-to-nearest mode can be easily implemented with minimal 
implementation cost (e.g., small amount of additional hardware and/or software) and 
marginally increases the processing time (if at all) for an arithmetic operation utilizing 

30 conventional flush-to-zero mode. In accordance with a specific implementation of the 
flush-to-nearest mode, a denormalized number is either flushed to zero or set to a 
particular set value (e.g., -a m i n , +a m i n ) depending on the value of the denormalized number 
and a threshold value. Other implementations and variants of the flush-to-nearest mode 
are described below. 
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An embodiment of the invention provides a floating-point unit (FPU) 



configurable to perform floating-point operations and having an operand processing 
section coupled to an operand flush section. For each floating-point operation, the 
operand processing section receives and processes one or more input operands to provide 
5 a preliminary result. The operand flush section determines whether the preliminary result 
falls within one of a number of ranges of values, and sets the preliminary result to one of 
a number of "set" values if the preliminary result falls within one of the ranges. The 
ranges typically cover denormalized numbers. 



10 preliminary result to one of two set values if it falls within one of two ranges of values. 
The first range can include values greater than zero and less than half of a positive 
minimum normalized number (i.e., 0 < y < +a m i n /2), and the second range can include 
values equal to or greater than half of the positive minimum normalized number and less 
than the positive minimum normalized number (i.e., +a m i n /2 < y < +a m i n ). The positive 

1 5 preliminary result can be set to zero if it falls within the first range and to +a m i n if it falls 
within the second range. Similar processing can be performed for a negative preliminary 
result such that it is flushed zero if it falls within the range (0 > y > -a m i n /2) and set to -a m i 
if it falls within the range (-a m i n /2 > y > -a m i n ). The ranges are typically selected for 
mathematical accuracy or for some other reasons. The ranges may also be selected such 

20 that the determination of whether the preliminary result falls within one of the ranges can 
be achieved by simply checking the exponent of the preliminary result. 



section and an exponent processing section. The mantissa processing section receives 
and processes one or more mantissas to provide a preliminary result mantissa, and the 

25 exponent processing section receives and processes one or more exponents to provide a 
preliminary result exponent. Depending on whether the result of the floating-point 
operation falls within the first or second range, the preliminary result mantissa can be 
respectively set to a first (e.g., 0.0. . .) or a second (e.g., 1.0. . .) mantissa value, and the 
preliminary result exponent can likewise be set to a first (e.g., 0) or a second (e.g., E m i n ) 

30 exponent value. Since the integer bit of the mantissa is hidden in some representations 
(such as for the IEEE standard), the action to set the mantissa to 0.0. . . or 1 .0. . . may be 
omitted. 



In a specific implementation, the operand flush section sets a positive 



The operand processing section typically includes a mantissa processing 



4 



Attorney Docket No.: 103-1 
Client l^fcence No.: 0099.00US 




Another embodiment of the invention provides a floating-point processor 



that includes a memory unit, an instruction dispatch unit, and a floating-point unit (FPU). 
The memory unit stores instructions and (possibly) operands for floating-point operations. 
The instruction dispatch unit retrieves instructions from the memory unit and provides the 
instructions to the FPU. In accordance with the received instructions, the FPU performs a 
floating-point operation on one or more input operands to provide a preliminary result, 
which is flushed to zero if it falls within a first range of values and set to a (plus or minus) 
set value if it falls within a second range of values. Again, the first range may be defined 
to include 0 < |y| < +a m i n /2, the second range may be defined to include +a min /2 < |y| < 
a m in, and the set value may be selected to be +a m i n or -a m j n depending on whether the 
preliminary result is positive or negative, respectively. Again, the ranges are typically 
selected for mathematical accuracy and may also be selected such that the determination 
on whether the preliminary result falls within the first or second range can be made by 
simply checking the exponent of the preliminary result. 



additional units to further process the operands. The FPU and floating-point processor 
can be incorporated within a microprocessor or other hardware structure, and can also be 
described and/or implemented using hardware design languages (e.g., Verilog). 



performing a floating-point operation. In accordance with the method, one or more input 
operands for the floating-point operation are received and processed to provide a 
preliminary result. A determination is then made whether the preliminary result falls 
within a first or second range of absolute values, which can be defined as described 
above. The preliminary result is flushed to a first value (e.g., zero) if falls is within the 
first range and set to a second value (e.g., +a m j n , +a min ) if it falls within the second range. 
Again, the determination on whether the preliminary result falls within the first or second 
range may be made by checking the exponent of the preliminary result. Various features 
described above also apply. 



are used to design and manufacture various embodiments of the invention. The invention 
further provides computer data signals embodied in a transmission medium and 
comprised various embodiments of the invention. 



The FPU and floating-point processor described above typically include 



Yet another embodiment of the invention provides a method for 



The invention also provides computer program products that implement or 



Attorney Docket No.: 103-1 
Client J^fcence No.: 0099.00US 




The foregoing, together with other aspects of this invention, will become 
more apparent when referring to the following specification, claims, and accompanying 
drawings. 



Fig. 1 is a block diagram of an embodiment of a processor that 
incorporates the operating mode of the invention; 

Fig. 2 is a block diagram of an embodiment of a portion of a floating-point 
unit (FPU) within the processor of Fig. 1; 



positive denormalized numbers; 

Fig. 5 is a diagram that graphically represents a reciprocal operation for 
some large numbers; 

Fig. 6 is a flow diagram of a specific implementation of the flush-to- 
nearest mode of the invention; 

Fig. 7 is a block diagram of an embodiment of a floating-point unit that 
can be used to implement the flush-to-nearest mode of the invention; and 

Fig. 8 is a logic diagram of a specific implementation of an exponent flush, 
unit and a mantissa flush unit shown in Fig. 7. 



incorporates the operating mode of the invention. Processor 100 includes a floating-point 
unit (FPU) 1 10, an instruction dispatch unit (IDU) 120, a memory 122, a load store unit 
(LSU) 130, and an integer execution unit (DCU) 140. FPU 110 implements a floating- 
point instruction set and performs floating-point computations, as directed by IDU 120. 
IDU 120 decodes a sequence of instructions, dispatches floating-point instructions to FPU 
110, and keeps track of the state of each dispatched floating-point instruction, resource 
and register dependencies, and the possibility of bypassing a resultant operand to the next 
FPU instruction. LSU 130 interfaces with other elements (i.e., internal or external to 
processor 100) and exchanges data with FPU 110. For example, operands are loaded 



BRIEF DESCRIPTION OF THE DRAWINGS 



Fig. 3 shows a representation of a normalized floating-point number; 
Fig. 4A is a diagram of a line that graphically represents all real numbers; 
Fig. 4B is a diagram of a line that graphically represents the range of 



DESCRIPTION OF THE SPECIFIC EMBODIMENTS 



Fig. 1 is a block diagram of an embodiment of a processor 100 that 
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from LSU 130 to FPU 110 and results are provided from FPU 1 10 to LSU 130. IXU 140 
performs integer computations, and is able to exchange data with FPU 110. 

In the embodiment shown in Fig. 1, a memory unit 122 couples to IDU 
120 and provides instructions (and possibly operands) for FPU 110. Memory unit 122 
5 can be a random access memory (RAM), a read only access memory (ROM), a Flash 
memory, or other type of memory. Memory unit 122 may be located within processor 
100 (as shown in Fig. 1) or external to processor 100. 

Fig. 1 also shows a block diagram of a specific embodiment of FPU 110, 
which includes a floating-point register file (FPR) 152, a floating-point pipe file (PIP) 
10 154, a floating-point multiplier (MUL) 156, a floating-point adder (ADD) 158, a floating- 
point exponent unit (EXP) 160, a floating-point divide/square-root unit (DSQ) 162, and a 
floating-point control unit (CTL) 164. The functions performed by these units are briefly 
j5j described below. 

FPR 152 includes a number of read ports and a number of write ports. The 

is? ; 

Si 15 read ports can be configured to read, for example, up to three operands for an arithmetic 
operation and one operand for a store operation. The write ports can be configured to 
write, for example, one operand for an arithmetic operation or a load operation. 

PIP 154 couples to FPR 152 and also interfaces with LSU 130 and IXU 
ill 140. For each instruction, PIP 154 selects and receives operands from FPR 152, a load 

20 pipe file, a result pipe file, or a ROM. PIP 154 then unpacks the received operands (e.g., 
l2| from an IEEE-compliant format) into an internal data format recognized and operable by 

the processing units within FPU 110. PIP 154 also packs the results from FPU 110 into a 
data format (e.g., IEEE-compliant format) required by the external circuitry and provides 
the packed results to FPR 152. 
25 MUL 156 couples to PIP 154 and executes floating-point multiply 

instructions as well as the multiply portion of compound instructions such as a multiply- 
add (MADD) instruction. MUL 156 receives the mantissas of two operands from PIP 
1 54. In an embodiment, MUL 1 56 is implemented with a multiplier array that generates a 
set of sum and carry outputs having full precision. The sum and carry are provided to, 
30 and combined in a carry-propagate adder (CPA) to generate a precise, unrounded 

resultant mantissa. The lower bits of the resultant mantissa are logically combined to 
form a "sticky" bit that, in combination with the round bit and the current rounding mode, 
are used to generate rounding information to be sent to a subsequent rounding unit within 
MUL 156. MUL 156 is described in further detail below. 



Is 
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ADD 158 couples to MUL 156 and PIP 154. ADD 158 executes floating- 



point add and subtract instructions as well as the add portion of compound instructions 
such as MADD. ADD 158 receives two operands and performs floating-point magnitude 
addition/subtraction using, for example, a prescale adder (PSA) operated in parallel with a 
5 massive cancellation adder (MCA). The final output is selected from one of the adders 
and provided to PIP 154 for storage. 



exponent difference of the subject operands and the first few bits of the mantissas of such 
operands (e.g., the integer bit(s) and the most significant bit of the fractional portion). In 

10 a specific implementation, the PSA is always used for addition and selectively used for 
subtraction when the result is guaranteed to be no less than 0.1000000 before rounding 
and normalization, and the MCA is selectively used for subtraction when the result is 
guaranteed to be less than 1.0 and may have a large number of leading zeros before 
rounding and normalization. Operands that satisfy both conditions may be processed by 

15 either the PSA or MCA. Often, this selection criteria results in selection of the PSA for 
additions and subtractions when the difference between operand exponents is greater than 
two, and selection of the MCA for subtractions when the difference between operand 
exponents is less than or equal to two. 



20 exponents of the operands from PIP 154, estimates an output exponent based on the 

operation being performed and the exponents of the operands, and adjusts the estimated 
output exponent based on the rounding and normalization performed on the mantissa 
(e.g., by MUL 156 and ADD 158). EXP 160 may also perform other functions such as 
overflow/underflow prediction and checks. 

25 In an embodiment, EXP 160 includes an exception prediction unit (EPU) 

170 that functions to generate an exception prediction signal when there is a possibility 
that an arithmetic computation performed by FPU 110 would cause an overflow or 
underflow exception. The exception prediction signal is provided to CTL 164 (not shown 
in Fig. 2), which reports the exception to IDU 120. EDU 120 then initiates or directs the 

30 required responsive processing. For example, upon receiving the exception prediction 

signal, IDU 120 can cause other processing units within processor 100 to temporarily halt 
operation. If it later turns out that the computation does not actually cause an overflow or 
underflow exception, IDU 120 is so notified, and the other processing units resume 



In an implementation, selection of which adder to use is based upon the 



EXP 160 couples to PIP 154 and ADD 158. EXP 160 receives the 
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operation. EPU 170 is described in detail in the aforementioned U.S. Patent Application 
Serial No. 09/363,638. 

In accordance with the invention, EXP 160 further includes an exponent 
flush unit (EFU) 180 that functions to set the exponent to one of a set of values, based on 
5 the value of the preliminary result from an arithmetic operation. EFU 180 implements the 
exponent portion of the flush-to-nearest mode of the invention and possibly other 
operating modes (e.g., flush-to-zero, IEEE-compliant modes, and so on). EFU 180 is 
described in further detail below. 

DSQ 162 couples to MUL 156 and operates on the mantissas for divide 
10 and square-root instructions. DSQ 162 is designed to implement a particular algorithm 
(e.g., a radix-2 SRT or a radix-4 SRT algorithm). 

CTL 164 is the control unit for FPU 110. CTL 164 receives floating-point 
□I instructions and directs the operation of various units within FPU 1 10 to generate the 
ft proper output. CTL 164 also receives the exception prediction signal from EPU 170 and, 
il 15 responsive to the signal, initiates or directs the required processing by other units. 
|| Fig. 2 is a block diagram of an embodiment of a floating-point unit (FPU) 

200 that is capable of performing add/subtract, multiply, and multiply-add (Madd) 
operations, and at least a portion (e.g., exponent computation) of divide, reciprocal, and 
square root operations. FPU 200 depicts a portion of FPU 110 shown in Fig. 1, and 
20 includes a multiplier unit 202 coupled to an adder unit 204. Support circuitry, such as 
that shown in Fig. 1, is not shown in Fig. 2 for simplicity. 

In the implementation shown in Fig. 2, multiplier unit 202 unit includes a 
multiplier array 210, a carry-propagation adder (CPA) 212, a rounding unit 216, and an 
exponent combination unit 220. Adder unit 204 includes the remaining units shown in 
25 Fig. 2. FPU 200 includes several features that increase accuracy, simplify the hardware 
design, and improve operational performance, as described below. 

At any given moment, FPU 200 can be configured to perform one of at 
least three different operations including addition/subtraction, multiplication, and Madd. 
These operations are expressed by the following: 

30 Fd = Fs±Ft , 

Fd - Fs^ Ft , and 
Fd = ±{{Fs*Ft)±Fr) , 



15! 
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where Fd is the resultant output and Fr, Fs, and Ft are three input operands. The Madd 
operation can be further decomposed into the following operations and their 
corresponding names: 

Fd = (Fs *Ft) + Fr , MADD 

5 Fd = (Fs • Ft)- Fr , MSUB 

Fd = -((Fs • Ft) + Fr) 9 N MADD 

Fd = -{{Fs • Ft) -Fr), N MSUB, 

As shown in Fig. 2, the mantissas, Mant S and Mant T, for two operands 
are provided to multiplier array 210. Multiplier array 210 can implement, for example, a 
10 Booth or modified Booth algorithm, and can include partial product generation logic and 
a number of carry-save adders. The partial product generation logic produces partial 
13 products based on the mantissas. The carry-save adders add a number of partial products 
m together and send the outputs to other carry-save adders in a tree-like fashion until only 

two numbers are left, the final sum and carry. In a specific implementation, the carry- 
15 save adders take in four terms and combine them into two, but other configurations are 

t = j 

rh possible. 

^ Multiplier array 210 thus multiplies the two operands and provides the 

Ml 

\!\ product in sum-and-carry format to CPA 212. CPA 212 combines the sum and carry and 

ill 

\T\ . provides the resultant mantissa to rounding unit 216 that processes the mantissa based on 
9f 20 a specified operating mode of FPU 200. The operation of rounding unit 216 is further 

described below. The processed mantissa comprises the output mantissa from multiplier 

unit 202. 

The exponents, Exp S and Exp T, of the two operands are provided to 
exponent combination unit 220 that combines the exponents for a multiply operation. 
25 The combined exponent from unit 220 comprises the output exponent from multiplier unit 
202. 

For improved performance (i.e., faster operating speed), adder unit 204 
includes a prescale adder (PSA) 240 and a massive cancellation adder (MCA) 242 
operated in parallel. Depending on the characteristics of the operands, the output from 
30 either PSA 240 or MCA 242 is selected. 

To perform a floating-point addition, the mantissas of the two operands are 
typically aligned by shifting one of the mantissa and adjusting its exponent until the 
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exponents of the two operands are equal. The mantissas are then combined (e.g., added 
or subtracted), and the resultant mantissa is normalized. The number of shifts prior to the 
combination can be large (e.g., when adding a large number with a small number), and 
the number of shifts after the combination can also be large (e.g., when subtracting two 
operands having similar magnitudes). PSA 240 and MCA 242 are designed to efficiently 
process most input conditions, including these two extreme conditions. 



232a and 232b. The mantissas for operands R and T are provided to a MUX 230 that 
selects one of the mantissas, based on the operation to be performed, and provides the 
selected mantissa to MUXes 232a and 232b. MUX 232a selects the mantissa of the 
smaller operand and provides the selected mantissa to a right-shift unit 252. MUX 232b 
selects the mantissa of the larger operand and provides the selected mantissa to a CPA 



provided to a MUX 222 that selects one of the exponents based on the operation to be 
performed. The selected exponent from MUX 222 and the combined exponent from unit 
220 are provided to an exponent calculation unit 224 that determines the difference 
between the two exponents, E d ifr, and a preliminary result exponent, Ep re . Determination 
of the preliminary result exponent is dependent upon the arithmetic equation being 
performed, and is further described in the aforementioned U.S. Patent Application Serial 
No. 09/363,638. The preliminary result exponent (e.g., the larger exponent when 
performing an add operation) is provided to an exponent adjustment unit 226 and the 
exponent difference is provided to right-shift unit 252 that shifts the mantissa of the 
smaller operand to the right by the indicated exponent difference. The shifted mantissa is 
provided to CPA 254 that combines the two mantissas and provides a combined mantissa 
to a rounding unit 256. Rounding unit 256 rounds the combined mantissa and provides 
the rounded mantissa to a normalization unit 258. 



xxxx) or a non-normalized format (e.g., lxx.xxx—xxxx, lx.xxx—xxxx, or 0.1xxx«xxxx). 
Normalization unit 258 normalizes the result to the Ol.xxx—xx format by performing a 1- 
bit or 2-bit right-shift or left-shift, if necessary. The exponent is adjusted by exponent 
adjustment unit 226, as necessary, based on the normalization performed by 
normalization unit 258. 



For PSA 240, the mantissa from rounding unit 216 is provided to MUXes 



254. 



The exponents of operands R and T (Exp R and Exp T, respectively) are 



The mantissa from CPA 254 can be in a normalized format (i.e., Ol.xxx— 
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MCA 242 includes a CPA 262, a leading zero anticipator (LZA) 264, and 
a left-shift unit 266. For MCA 242, the mantissas from rounding unit 216 and MUX 230 
are provided to MUXes 232c and 232d. MUXes 232 facilitate a small shift of one of the 
mantissas, based on the exponent difference, to align the mantissas. MUXes 232 are also 
5 used to select one of the mantissas for inversion in subtraction operations (the inverter is 
not shown in Fig. 2 for simplicity). The outputs from MUXes 232c and 232d are 
provided to CPA 262 and LZA 264. CPA 262 combines the two mantissas and provides 
a combined mantissa to left-shift unit 266. LZA 264 anticipates the number of leading 
zeros in the resultant mantissa, based on the input operands. The output from LZA 264 is 
10 encoded into a control signal indicative of the number of shifts for left-shift unit 266. The 
control signal is also provided to exponent adjustment unit 226 to adjust the exponent. 

The outputs from normalization unit 258 and left-shift unit 266 are 
provided to a MUX 246 that selects the output from PSA 240 or MCA 242 as the output 
mantissa from FPU 200. The adjusted exponent from unit 226 comprises the output 
1 5 exponent from FPU 200. The operation of FPU 200 is further described in the 
aforementioned U.S. Patent Application Serial No. 09/364,514. 

Fig. 3 shows a representation of a normalized floating-point number. The 
representation includes a sign bit 310, a mantissa (mant) 312, and an exponent (exp) 314. 
A floating-point number (y) can be expressed as: 

20 y = (~l) sign • mant • 2 exp . Eq. (1) 

A normalized floating-point number is represented by a mantissa having a 
one (1) to the left of the binary point (i.e., the integer portion of the mantissa) and a 
1 .xxx — xx format, where each "x" represents one bit that is either a one or a zero. As 
, defined by the IEEE standard, the fractional portion "xxx~xx" represents 23 bits after the 
25 binary point for a normalized single-precision number and 52 bits after the binary point 
for a normalized double-precision number. The normalized mantissa has a range of 
between 1.0 and 2.0 (i.e., 1.0 < mantissa < 2.0). 

The IEEE standard defines the representation for floating-point numbers. 
For normalized numbers, the IEEE standard mandates storing only the fractional portion 
30 of the mantissa (i.e., the "xxx--xx" portion in Fig. 3). The leading one (1) to the left of 
the binary point is implicit and not stored. 

Fig. 4A is a diagram of a line that graphically represents all real numbers. 
It should be noted that this line representation is not drawn to scale. Generally, real 
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numbers range from negative infinity (-oc) to positive infinity (+oc). In the line 
representation shown in Fig. 4A and for the floating-point expression shown in equation 
(1), numbers greater than 0.0 (i.e., y > 0.0) are represented by a positive sign bit (i.e., sign 
= 0) and numbers less than 0.0 (i.e., y < 0.0) are represented by a negative sign bit (i.e., 
5 sign = 1). Numbers having an absolute value greater than 1.0 (i.e., |y| > 1.0) are 

represented by a positive exponent (i.e., exp > 0) and numbers having an absolute value 
less than 1 .0 (i.e., |y| < 1 .0) are represented by a negative exponent (i.e., exp < 0). 

For floating-point representations having finite resolution (e.g., 32-bit 
single-precision and 64-bit double-precision floating-point numbers), only numbers 
10 having absolute values within a particular range of values can be represented as 

normalized numbers using the expression shown in equation (1). This range is defined by 
a positive maximum normalized value, +a max , and a positive minimum normalized value, 
+a m i n . For single-precision numbers, +amax is slightly smaller than 2 +l28 and +a min = 2" 126 . 
Some numbers having absolute values between zero and the positive minimum 

j 15 normalized value (i.e., 0.0 < y < +a min ) can be represented as denormalized numbers, 

which are further described below. Special representations are used for zero (0.0) and for 
some numbers greater than the maximum normalized value (e.g., positive infinity and 
negative infinity). 

Table 1 lists the representations for some of the normalized and 

i 20 denormalized numbers. As suggested in Table 1, successively smaller normalized 

numbers are obtained by decrementing the mantissa from 1.111—11 down to 1 . 000—00. 
The next successively smaller normalized number is obtained by resetting the mantissa to 
1.111 — 11 and decrementing the exponent by one. This continues until the mantissa is 
1.000-00 and the exponent is at the minimum value, Emin, used for representing 

25 normalized numbers. The maximum denormalized number is defined by a mantissa of 
0.1 1 1—1 1 and an exponent of Emm- Similarly, successively smaller denormalized 
numbers are obtained by decrementing the mantissa from 0.1 1 1 — 1 1 down to 0. 000—01. 

Table 1 - Floating Point Number Representation 



Number 


Floating Point 
Representation 


± infinity 


±l.lll-ll*2 Emax+1 


± max normalized number (± amax) 


± 1.111 — ll*2 Emax 
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± min normalized number (± a m j n ) 


± 1.000~00»2 Emin 


± max denormalized number 


± 0.111— ll»2 Emin 






± min denormalized number 


± 0.000-0 l»2 Emin 


± zero (± 0) 


± 0.000--00»2° 



As noted above, the IEEE standard defines a number of "rounding" modes 
for converting a preliminary result from an arithmetic operation (which is usually 
represented in an internal data format) into an IEEE-compliant final result. The rounding 
modes define by the IEEE standard include the "round-to-nearest", "round-to-zero", 
"round-to-positive-infinity", and "round-to-negative-infinity" modes. 

In the round-to-nearest mode, the preliminary result is calculated to be the 
closest representable floating-point number, which can be either bigger or smaller than 
the preliminary result. In the round-to-zero mode, the preliminary result is calculated to 
be the nearest representable floating-point number toward zero. Thus, the preliminary 
result is calculated to be the next smaller floating-point number if it is a positive number, 
and calculated to be the next larger floating-point number if it is a negative number. In 
the round-to-positive-infinity mode, the preliminary result is calculated to be the next 
larger floating-point number regardless of whether it is a positive or negative number. 
And in the round-to-negative-infinity mode, the preliminary result is calculated to be the 
next smaller floating-point number, again regardless of whether it is a positive or negative 
number. 

The IEEE standard defines performing the rounding at a particular bit 
position of a floating-point number. This requirement ensures a consistent result from 
different implementations of a floating-point processor. To obtain an IEEE-compliant 
result, the preliminary result is normalized and then rounded. 

Many conventional floating-point processors also support a "flush-to-zero" 
operating mode in which denormalized preliminary results are flushed to zero. The 
processing for a flush-to-zero operation is typically simpler to implement than the 
processing for a rounding operating and can also be executed in less time. However, the 
flush-to-zero mode typically has accuracy loss as a result of the flush. Many embedded 
processors implement the flush-to-zero mode (and not the rounding modes described 
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above) because the loss in accuracy can be tolerated and the reduced complexity and 
improved operating speed are highly desirable. 

The invention provides an operating mode that has improved accuracy 



over the flush-to-zero mode. In one embodiment, the operating mode of the invention can 
5 be easily implemented with minimal implementation cost (e.g., small amount of 

additional hardware and/or software) and marginally increases the processing time (if at 
all) for an arithmetic operation utilizing conventional flush-to-zero mode. Li accordance 
with this operating mode, which is also referred to as a "flush-to-nearest" mode, a 
denormalized number is either flushed to zero or set to a particular "set" value depending 

10 on the value of the denormalized number and a threshold value. The flush-to-nearest 

mode of the invention can provide more accurate final results for at least some arithmetic 
operations, as described below. In a preferred embodiment, the "set" values are 
predetermined and maintained in a conventional format (e.g., hard wired as a constant 
value, stored in a register as a changeable value, etc.) to avoid repeated calculations. 

15 In a specific implementation of the flush-to-nearest mode, the threshold 

value is selected to be the mid-point between zero and the minimum normalized number 
(i.e., threshold = ±a m i n /2) and the set value is selected to be the minimum normalized 
number (i.e., set = ±a m i n ). In this implementation, a positive denormalized number is 
flushed to zero if it is greater than zero and less than +a m i n /2 (i.e., 0 < y < +amin/2), and set 

20 to -t-a m in if it is equal to or greater than +a m m/2 and less than +a m i n (i.e., +a m i n /2 < y < 

+a m in). Similarly, a negative denormalized number is flushed to zero if it is less than zero 
and greater than -a m i n /2 (i.e., 0 > y > -a m i n /2), and set to -a m i n if it is equal to or less than 
-a m i n /2 and greater than -a mm (i.e., -a m i n /2 > y > -a min ). 



25 values (e.g., ±a m i n /4, ±a m i n /8, and others) and other set values (e.g., ±a m i n /4, ±3a m i n /4, and 
others) may be used. These other threshold and set values may provide better results for 
some applications. In these implementations, a denormalized number is flushed to zero 
if it falls within a first range of values and set to the set value if it falls within a second 
range of values. For example, a positive denormalized number may be flushed to zero if 

30 it is greater than zero and less than +a m j n /4 (i.e., 0 < y < +a min /4), and set to +a m i n if it is 
equal to or greater than +a m i n /4 and less than +a m i n (i.e., +a m j n /4 < y < +a m i n ). 



positive denormalized numbers, which is flanked on the left side by positive zero (+0.0) 



In other implementations of the flush-to-nearest mode, other threshold 



Fig. 4B is a diagram of a line that graphically represents the range of 
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and on the right side by the positive minimum normalized number, +a m j n , which is also 
referred to as positive min-norm. The mid-point between zero and +a m i n is shown by a 
hashed mark between these end points. As shown in Fig. 4B, a denormalized number 
greater than zero and less than +a m i n /2 is flushed to zero and a denormalized number equal 
to or greater than +amin/2 and less than +a m i n is set to +a m i n . 

The FPU can be designed to support a combination of operating modes. 
Table 2 lists some of the operating modes and their attributes. The FPU can be selected 
to operate in accordance with one of the supported operating modes, which can be 
selected based on, for example, a status register such as that described in the 
aforementioned U.S. Patent Application Serial No. 09/364,786. 



Table 2 



Operating Mode 


Implementation 
Complexity 


Acccuracy 


IEEE 


moderate 


good 


Flush-to-Zero 


simple 


poor 


Flush-to-Nearest 


simple 


improved over the 
Flushed-to-Zero 



Fig. 5 is a diagram that graphically represents a reciprocal operation for 
some large numbers. In a reciprocal operation, a number greater than 1 .0 is transformed 
into a number smaller than 1 .0. In accordance with the IEEE standard and for a single- 
precision floating-point number, an exponent range of 0 through +127 is used to represent 
numbers greater than 1.0 and an exponent range of -1 through -126 is used to represent 
numbers less than 1.0. The range of representable normalized numbers greater than 1.0 is 
thus greater than the range of normalized numbers less than 1.0. Consequently, some 
large normalized numbers near positive infinity will map to denormalized numbers in a 
reciprocal operation. 

For a reciprocal operation implemented using a Newton-Raphson 
algorithm, a number of iteration is performed to achieve the final result. For some large 
normalized numbers, the results of the first iteration will map to the denormalized number 
range. When operating in the conventional flush-to-zero mode, these intermediate results 
from the first iteration are flushed to zero, as indicated by the dashed line, and accuracy is 
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lost. However, when operating in the flush-to-nearest mode of the invention, 
intermediate results equal to or greater than +a m i n /2 and less than +a m i n will be set to +a m i n , 
and intermediate results equal to or less than -a min /2 and greater than -a m i n will be set to 
-a m i n . Subsequent iterations may then be able to provide the reciprocal results (e.g., a 
result of x as shown in Fig. 5). 

Thus, the flush-to-nearest mode of the invention increases the accuracy of 
at least some arithmetic operations (e.g., reciprocal, multiply, divide, and so on) in which 
the results of the operations may be denormalized numbers. 

For the specific implementation of the flush-to-nearest mode described 
above in which the threshold value is ±a m i n /2 and the set value is ±a m in, accuracy increases 
by 50% on the average over the conventional flush-to-zero mode. 

The threshold value and set value are typically selected to provide 
mathematical accuracy, and may further be selected such that the flush-to-nearest mode 
of the invention can be implemented with small implementation cost and minor (if any) 
degradation in operating speed. For example, if the threshold value is selected to be 
related to the minimum normalized number by a power of two (e.g., ±a m j n /2), the 
determination for whether a preliminary result is equal to or greater than the threshold 
value can be easily achieved. Referring back to Table 1, the positive minimum 
normalized number is expressed as: 

+a min = 1.000--00«2 Emin . 

A positive threshold value of +a m j n /2 can thus be expressed as: 

+a min /2 = 1.000~00#2 Emin l . 

From the above expressions, it can be observed that a positive preliminary result having a 
value equal to or greater than +a m j n /2 and less than +a m i n (i.e., +ami n /2 < y < +amj n ) will 
have an exponent value of E m i n - 1 . Thus, to determine whether the preliminary result is 
equal to or greater than +a m i n /2 and less than +a m i m it is only necessary to check the 
exponent of the preliminary result to determine whether it is equal to E m i n - 1 . There is no 
need to check the mantissa of the preliminary result. A negative preliminary result differs 
from a positive preliminary result only in the sign bit, and can be checked in similar 
manner as for a positive preliminary result. 

From the above description, it can be observed that the selection of a threshold value of 
±a m in/2 corresponds to a detection of E m in - 1 in the preliminary result exponent; Emi n - 1 
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is thus an exponent comparison value to be used to determine whether a positive 
preliminary result is equal to or greater than +a m i n /2 and less than +a m i n . A different 
threshold value may correspond to a detection of a different exponent value. For 
example, selection of a threshold value of +a m i n /4 corresponds to a detection of E m j n - 2 
5 (and possibly Emi n - 1) in the exponent. By selecting the threshold value to be a power of 
two smaller than a m i n (e.g., +a m i n /2, +a m i n /4, and so on) the detection process for one end 
of the range (e.g., +a m i n /2 < y) is greatly simplified since only the exponent needs to be 
checked and the mantissa can be ignored. However, other threshold values not related to 
+a m in by a power of twos can also be used and are within the scope of the present 
10 invention. Proper selection of the equality (i.e., +a m i n /2 < y, and not +a m i n /2 < y) also 
greatly simplifies the detection process. 

Thus, by selecting the threshold value to be +a m i n /2 and the upper end of 
the range to be +a m in, where the threshold value and the upper end is related by a factor of 
two, and by properly choosing the equality conditions (e.g., +a m in/2 < y < +a m i n ), the 
1 5 detection of whether a positive preliminary result falls within the range can be made by 
simply comparing the exponent of the positive preliminary result against one exponent 
comparison value (e.g., E m j n - 1). These design conditions greatly simplify the 

!L implementation of the flush-to-nearest mode of the invention. However, these design 

RJ 

111 conditions are not required by the invention. 

Ill 

i s j 20 The specific implementation of the flush-to-nearest mode described above 

Nj uses one threshold value (e.g., ±a m i n /2) and one set value (e.g., ±a m i n ). Generally, the 

RS 

invention can be extended to include any number of threshold values and any number of 
set values. Specifically, N threshold values in combination with zero (or some other 
value) can be used to define N+l ranges of values. The preliminary result can be flushed 

25 to zero or set to one of N set values, depending on whether the preliminary result falls 
into one of the N+l ranges. If zero is considered one of the set values, the result is thus 
effectively be set to one of N+l set values, depending on which one of the N+l ranges it 
falls into. For example, a flush-to-nearest quartile mode can be implemented in which a 
preliminary result is set to the nearest quartile (e.g., 0, ±a m m/4, ±a m i n /2, ±3a m i n /4, or ±a m i n ). 

30 The implementation for the flush-to-nearest quartile will be different than that described 
above for the flush-to-nearest 0 or ±a m i n , and may entail checking the exponent as well as 
the mantissa of the preliminary result (e.g., for ±3a m i n /4). Other variations of the flush-to- 
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nearest mode of the invention can be contemplated and are within the scope of the present 
invention. 

Fig. 6 is a flow diagram of a specific implementation of the flush-to- 
nearest mode of the invention. This flow diagram is for the specific implementation 
5 described above in which a preliminary result is flushed to zero if its absolute value is 
within the range of 0 < |y| < +a min /2, and set of +a min or -a min , depending on the sign bit, if 
its absolute value is within the range of +a m i n /2 < |y| < +a m i n . Initially, at step 610, the 
exponent(s) for a particular arithmetic operation are processed in the normal manner to 
determine a preliminary result. A test is then performed to detect underflow and to 
10 determine whether the exponent of the preliminary result is equal to Emin - 1, at step 612. 
Underflow can be predicted based solely on the one or more exponents of the one or more 
operands for the arithmetic operation being processed. For example, underflow can be 
\2\ predicted based on one exponent for a square root or reciprocal operation, two exponents 
sl{ for an add, subtract, multiply, or divide operation, and three exponents for a multiply-add 

y\ 15 operation. Underflow prediction is described in detail in the aforementioned U.S. Patent 
m Application Serial No. 09/363,638. 

W Based on the result of the underflow test and the exponent check 

% performed in step 612, a set of determinations is made. If an underflow is detected and 

{?! the exponent is equal to E min - 1, at step 620, the preliminary result is set to ±a m i n5 at step 

j Jj 20 622. If an underflow is detected and the exponent is not equal to Emin - 1 , at step 630, 
13 indicating a preliminary result having an absolute value that is less than +a m i n /2, the 
^ preliminary result is flushed to zero, at step 632. Otherwise, if no underflow is detected, 

at step 640, the normal result is provided, at step 642. 

of the steps shown in Fig. 6 are normally performed in floating- 
25 point processors, such as?he^e described in the aforementioned Patent Application Serial 
Nos. 09/363,638, 09/364,514, 09/5^4^12, 09/334,927, 09/364,786, 09/364,789, and 
09/364,787. To implement the flush-to-neafrs$Jmode of the invention, step 612 is 
modified to include the check to determine whetheHl^e exponent is equal to E m i n -1 and 
steps 620 and 622 are added. It can be noted that in these^tegs, only the preliminary 
30 exponent, Exp, is checked to determine whether it is equal to EminVJ. The additional 
steps can be easily achieved with small implementation cost and can be^s^ecuted with 
small or no additional processing delay. 
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Fig. 7 is a block diagram of an embodiment of a floating-point unit (FPU) 



700 that can be used to implement the flush-to-nearest mode of the invention. FPU 700 
includes an exponent processing section that processes the exponents of the input 
operands and a mantissa processing section that processes the mantissas of the input 
5 operands. In the embodiment shown in Fig. 7, the exponent processing section includes 
an exponent processing unit 710, an underflow/overflow prediction unit 712, an exponent 
adjustment unit 714, and an exponent flush unit 716. In this embodiment, the mantissa 
processing section includes a mantissa processing unit 720 and a mantissa flush unit 722. 
For simplicity, FPU 700 is shown processing two operands (e.g., for an add, subtract, 
10 multiply, or divide operation) but may be configured to process fewer operands (e.g., for 
a square root or reciprocal operation) or more operands (e.g., for a multiply-add 
operation). 



coupled to an operand flush section. The operand processing section includes exponent 

15 processing unit 710, underflow/overflow prediction unit 712, exponent adjustment unit 
714, and mantissa processing unit 720. The operand flush section includes exponent 
flush unit 716 and mantissa flush unit 722. Other embodiments of an FPU that can 
implement the flush-to-nearest mode of the invention can be designed and are within the 
scope of the present invention. 

20 As shown in Fig. 7, the exponents of the two operands, Exp A and Exp B, 

are provided to exponent processing unit 710 that processes the received exponents to 
generate a set of intermediate exponent values. Exponent processing unit 710 generally 
includes exponent combination unit 220 and exponent calculation unit 224 in Fig. 2. The 
intermediate exponent values may include, for example, a preliminary result exponent 

25 (Ep re ), an exponent difference (Edifr) between the received exponents, an overflow 

threshold value (OTV), and an underflow threshold value (UTV). The preliminary result 
exponent, Ep re , is determined based on the type of operation being performed and the 
values of the received exponents, Exp A and Exp B. The preliminary result exponent, 
Ep re , and the exponent difference, E d iff, are used for the mantissa processing and for the 

30 underflow and overflow exception prediction. The computations to derive the 

intermediate exponent values are described in further detail in the aforementioned U.S. 
patent application Serial No. 09/363,638. 



prediction unit 712 that, among other functions, predicts whether an underflow is likely to 



FPU 700 can also be viewed as including an operand processing section 



The intermediate exponent values are provided to underflow/overflow 



20 



Attorney Docket No.: 103-1 
Client^^rence No.: 0099.00US 

occur. The preliminary result exponent, Ep re , along with a signal from mantissa 
processing unit 720 are provided to exponent adjustment unit 714 that adjust the exponent 
(e.g., by -2, -1, 0, or +1) based on the result of the mantissa processing. Exponent 
adjustment unit 714 generally corresponds to exponent adjustment unit 226 in Fig. 2. 
5 The preliminary result exponent, Exp D, from exponent adjustment unit 

714 and a "constant" exponent, Exp E, are provided to exponent flush unit 716. If Exp D 
is equal to Emm - 1, indicating that the preliminary result falls within the range +a m i n /2 < 
|y| < +a m i m then the output exponent, Exp F, is set to Emm to set the final result to ±a m i n . 
Otherwise, if Exp D is less than Emj n - 1, indicating that the preliminary result falls within 
10 the range 0 < |y| < +a m i n /2, then the final exponent, Exp F, is set to E m i n - 1 to flush the 
output result to zero. Exp E thus has a value of either E m i n or E min - 1 depending on 
whether the result is to be set to ±a m i n or flushed to zero. Exponent flush unit 716 may 
3 further provide a status signal indicative of whether the preliminary result exponent is 
|| being flush to zero or set to Emi n . The status signal is provided to, and used by, mantissa 

a. ! 

15 flush unit 722 to flush or set the preliminary result mantissa, 
l! Mantissa processing unit 720 receives and processes the mantissas of the 

li two operands, Mant A and Mant B, to generate a preliminary result mantissa, Mant D. 
Lf* Mantissa processing unit 720 generally includes the units to the right of line 206 in Fig. 2. 
J! The preliminary result mantissa, Mant D, is then provided to a mantissa flush unit 722 
jj 20 that operates on similar principle as exponent flush unit 716. Specifically, the 
g preliminary result mantissa is flushed to a value corresponding to zero if the preliminary 
result falls within the range 0 < |y| < +ami n /2, and set to a value corresponding to ±a min if 
the preliminary result falls within the range +a m i n /2 < |y| < +amj n . The determination to 
flush the mantissa or to set it to a value corresponding to ±a m i n can be achieved using the 
25 status signal provided by exponent flush unit 716. 

In FPU 700, underflow/overflow prediction unit 712 is used to provide 
improved performance (i.e., reduced processing delay). However, the flush-to-nearest 
mode of the invention can be implemented without underflow/overflow prediction unit 
712, and this is within the scope of the present invention. 
30 Fig. 8 is a logic diagram of a specific implementation of exponent flush 

unit 716 and mantissa flush unit 722 in Fig. 7. A Flush-to-Zero signal indicative of the 
FPU operating in the flush-to-zero mode and a Flush-to-Nearest signal indicative of the 
FPU operating in the flush-to-nearest mode are provided to the input of an OR gate 812. 
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The output of OR gate 812, which is logic high if the FPU is operating in either the flush- 
to-zero or flush-to-nearest mode, is provided to an input of an NAND gate 814. An 
Underflow signal from underflow/overflow prediction unit 712 is provided to the other 
input of NAND gate 814. The output of NAND gate 814 comprises a FLUSHB signal 
5 and is logic low if an underflow is detected and the FPU is operating in either the flush- 
to-zero or flush-to-nearest mode. 

The FLUSHB signal is provided to one input of a number of flush circuits 
820a through 820n. Each of flush circuits 820a through 820n also receives a respective 
bit from the exponent or the mantissa. Within flush circuit 820a, the FLUSHB signal and 

10 the received bit from the exponent or mantissa are provided to the inputs of a NAND gate 
822a. The output of NAND gate 822a couples to an input of an inverter 824a, and the 
output of inverter 824a comprises the output bit for the exponent or mantissa. If the 
FLUSHB signal is logic low, indicating the need to flush the result, the output of NAND 
gate 822a is logic high and the output of inverter 824a is logic low. Each of flush circuits 

15 820b through 820n are similarly configured and operated. 

- A signal indicative of the condition (Exp = Emi n - 1 ) and the Flush-to- 
Nearest signal are provided to the input of an AND gate 828. The output of AND gate 
828 comprises a SET signal, which is logic high when the preliminary result is to be set to 
+a min or -a min . The FLUSHB signal, the SET signal, and the least significant bit (LSB) of 

20 the exponent are provided to a flush/set circuit 830. 

Within flush/set circuit 830, the FLUSHB signal and the exponent LSB are 
provided to the input of a NAND gate 832, similar to the other flush circuits 820. The 
output from NAND gate 832 is logic high if the exponent is to be flushed to zero. 

The FLUSHB signal and the SET signal are also respectively provided to 

25 the inverting and non-inverting input of a NAND gate 834. The output from NAND gate 
834 is logic low if both input conditions are true (i.e., the FLUSHB signal is logic low 
and the SET signal is logic high) and the exponent is to be set to E min . The outputs from 
NAND gates 832 and 834 are provided to the inputs of a NAND gate 836. If the output 
from NAND gate 834 is logic low, the output from NAND gate 836 is forced to logic 

30 high and the exponent is set to E min . Otherwise, the exponent LSB is flushed to zero if the 
FLUSHB signal is logic low, similar to the other exponent and mantissa bits. In the 
embodiment shown in FIG. 8, there is no need to flush the MSB of the mantissa to 0 (to 
flush the operand to ±0.0 or ±0.0. . .x 2°) or to set the MSB of the mantissa to 1 (to set the 
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operand to ±a m i n or ±1 .0. . .x 2 min ) since the integer bit in the mantissa is hidden in the 
IEEE standard. 

As can be seen for the specific implementation shown in Fig. 8, minimum 
additional circuitry is required to set the preliminary result to +a m j n or -a m j n . 

In accordance with the IEEE standard, a single-precision floating-point 
number is represented by a 32-bit value comprised of a sign bit, an 8-bit exponent, and a 
23-bit mantissa, and a double-precision floating-point number is represented by a 64-bit 
value comprised of a sign bit, an 1 1-bit exponent, and a 52-bit mantissa. Representations 
for normalized and denormalized, single and double-precision floating-point numbers are 
described in further detail in the aforementioned IEEE standard and in U.S. Patent 
Application Serial No. 09/364,5 12. 

In many floating-point processors, EEEE-compliant floating-point numbers 
are converted to, and processed using, an internal representation having a greater range 
such that large intermediate and preliminary results can be accurately represented. For 
example, a processor can be designed capable of internally representing and processing 
12-bit exponents that can accurately represent the 8-bit single-precision and the 1 1-bit 
double-precision exponents. 

Table 3 lists the true exponent values, their external representation in 
accordance with the IEEE standard for single-precision numbers, and their internal 
representation in accordance with a specific implementation using 12-bit internal 
exponent representation. In this implementation, the true exponent values are represented 
internally by inverting the most significant bit (MSB) of the external representation and 
sign extending the MSB. For example, for a true exponent value of zero, the single- 
precision IEEE representation is x7F (hex) or 0111,1111 (binary) and the internal 
representation is xFFF (hex) or 1111,1111,1111 (binary). Similarly, for a true exponent 
value of E min , the single-precision IEEE representation is xOl (hex) or 0000,0001 (binary) 
and the internal representation is xF81 (hex) or 1 1 1 1,1000,0001 (binary). 



Table 3 



Exponent Value 


Single-Precision 

' External 
Representation 


Single-Precision 

Internal 
Representation 


infinity 


xFF 


x07F 


Emax 


xFE 


x07E 
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0 


x7F 


xFFF 








Emin 


xOl 


xF81 


zero/denorm 


xOO 


xF80 



For the specific internal representation described above and shown in 
Table 3, E m i n for single-precision floating-point numbers is represented as xF81 (hex) and 
E min for double-precision floating-point numbers is represented as xCOl (hex). For this 
specific implementation, a denormalized number can be determined to be equal to or 
greater than +a min /2 and less than +a m j n if the internal exponent is xF80 for a single- 
precision floating-point number and xCOO for a double-precision floating-point number. 
For single-precision floating-point numbers, an underflow is detected if the exponent is 
xF80 or less. A check to determine whether the exponent is Emin - 1 can be performed by 
checking whether the upper five bits are all ones for single-precision and whether the 
upper two bits are all ones for double-precision. Thus, for the specific implementation 
described above and in conjunction with the underflow detection, not all bits in the 
exponent need to be check to determine whether the exponent is equal to Emin - 1 . 

For clarity, the invention has been described in the context of single- 
precision and double-precision floating-point representations that conform to the IEEE 
formats. However, the invention can be adopted for used with other floating-point 
representations, and this is within the scope of the invention. 

The FPUs described above can be implemented in various manners. For 
example, the FPUs can be implemented in hardware within a digital signal processor, an 
application specific integrated circuit (ASIC), a microprocessor, an embedded processor, 
and other hardware structures. 

In addition to implementations of the invention using hardware, the 
invention can also be embodied in a computer usable medium configured to store a 
computer-readable program code. The program code causes the enablement of the 
functions or fabrication, or both, of the invention disclosed in this specification. For 
example, this might be accomplished through the use of general programming languages 
(e.g., C, C++, and so on); hardware description language (HDL), including Verilog HDL, 
VHDL, AHDL (Altera hardware description language); or other programming and/or 
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circuit (i.e., schematic) capture tools available in the art. As a specific example, the 
Verilog simulator "VCS v.4.1.1" by Synopsys, Inc. may be used to simulate the 
invention. A book entitled "A Verilog HDL Primer" by J. Bhasker, Star Galaxy Pr., 1997 
provides greater detail on Verilog HDL, and is incorporated herein by reference in its 
entirety for all purposes. Using Verilog HDL, for example, a design may be described at 
a variety of levels, including switch-level, gate-level, register-transfer-level (RTL), and 
algorithmic level. In the program code implementation, Figs. 2, 6, and 7 can serve as an 
embodiment of a flow diagram. 



containing a computer usable (e.g., readable) medium including semiconductor memory, 
magnetic disk, optical disc (e.g., CD-ROM, DVD-ROM, and so on), and as a computer 
data signal embodied in a computer usable (e.g., readable) transmission medium (e.g., 
carrier wave or any other medium including digital, optical, or analog-based medium). 
As such, the code can be transmitted over communication networks, including the 
Internet and intranets. 



described above can be represented in a core that can be utilized in programming code 
and transformed to hardware as part of the production of integrated circuits. Also, other 
implementations of the invention (e.g., the FPUs described above) using a combination of 
hardware and software are also possible. Therefore, the embodiments expressed above 
are within the scope of the invention and should also be considered protected by this 
patent. 



enable any person skilled in the art to make or use the present invention. Various 
modifications to these embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other embodiments without 
the use of the inventive faculty. Thus, the present invention is not intended to be limited 
to the embodiments shown herein but is to be accorded the widest scope consistent with 
the principles and novel features disclosed herein, and as defined by the following claims. 



The program code can be disposed in any known article of manufacture 



It is understood that the functions accomplished by the invention as 



The foregoing description of the specific embodiments is provided to 
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